ASP Conference Series, Vol. **VOLUME***, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



On the Nature of X-Ray Flashes 



R. Mochkovitch, F. Daigne 

Institut d'Astrophysique de Paris, 75014 Paris, France 
C. Barraud, J.L. Atteia 

Observatoire Midi-Pyrenees, 31400 Toulouse, France 

Abstract. We have developed a toy model for internal shocks which has 
been used to generate a large number of synthetic GRBs in order to find 
in the parameter space the conditions which can lead to the formation of 
X-ray flashes. The key condition appears to be a small contrast of the 
Lorentz factor in the relativistic wind emitted by the central engine. 



1. Introduction 

The recently identified X-ray flashes (Kippen et al, 2002) have non thermal 
spectra with E p < 50 keV and weak gamma-ray fluxes (< 0.2 ph.cm _2 .s _1 in 
the 50 - 300 keV energy range). They were discovered by Beppo-SAX and are 
currently observed by HETE 2. They can be normal GRBs at very high red- 
shifts (which however does not seem compatible with their duration distribution 
similar to that of long bursts) or viewed off-axis (Yamasaki et al, 2002) but their 
soft spectra can also have an intrinsic origin. We have tested this possibility in 
the context of the internal shock model where it is supposed that a central en- 
gine generates a relativistic wind with a non uniform distribution of the Lorentz 
factor. 



2. A toy model for internal shocks 

The evolution of this relativistic wind should be followed with a fully hydro- 
dynamical calculation. This can be done (Daigne and Mochkovitch, 2000) but 
requires large amounts of computing time which prevents to consider a large 
number of cases and fully explore the parameter space. These detailed cal- 
culations have however shown that a simplified approach where the wind is 
represented by many shells which interact by direct collisions only (all pressure 
waves being neglected) can also produce satisfactory results (Kobayashi, Piran 
and Sari, 1997; Daigne and Mochkovitch, 1998). Going a step further we have 
developed for this study a toy model where internal shocks are limited to the 
collision of two shells of equal mass m. Shell 2 (Lorentz factor T2) is generated 
a time r after shell 1 (Lorentz factor Y\ < I^). The average power injected into 
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the wind in this two shell approximation is given by 



2 

e = — (ri + r 2 ) = mTc 2 (i) 



where M = 2m /t and T = ^{Y\ + Y-2) are the average mass loss rate and Lorentz 

p2p2 

factor. Shell 2 will catch up with shell 1 at the shock radius r s = 2cT r2 1 _3, i . 

1 2 1 1 

The two shells merge and the energy dissipated in the collision -Ediss is radiated 
with a characteristic broken power law spectrum. If the synchrotron process is 
responsible for the emission, the peak energy (maximum of uF u ) is 

E p ~E syn <xr s Br 2 e <xT s p x s ey (2) 

In the last term p s and e s are respectively the typical density and dissipated 
energy (per unit mass) in the comoving frame of the shocked material. The 
standard equipartition assumptions correspond to x = 1/2 and y = 5/2 but we 
consider below the possibility that x and y can have different values. This may 
be for example the case if the equipartition parameters «b an d a e which fix the 
fraction of dissipated energy injected into the magnetic field and the relativistic 
electrons are not constant but depend on p s or /and e s . 

The physical properties of the shocked material T s , p s and e s can be directly 
related to the wind parameters r, E, T and k = ^/Ti so that eq.(2) can be 
expressed as 

E P ^ r 2xf6x-l ^ 

where <p xy (K) is an increasing function of n for all reasonable values of x and y. In 
spite of the simplicity of the two shell approximation eq.(3) predicts a duration- 
hardness relation and a hardness-intensity correlation (HIC) as observed in real 
bursts. Another interesting (and surprising) consequence of eq.(3) is that E p is 
a decreasing function of F as long as x > 1/6. . Soft bursts are obtained with 
"clean fireballs" (large T) while winds with larger baryon loads ("dirty fireballs") 
lead to a harder spectrum because internal shocks occur closer to the source at 
a higher density. 

From the values of Efo ss , r, -EWeak and assuming a Band spectrum with 
low and high energy indices a = — 1 and (3 = —2.5 it is possible to estimate 
(for a given redshift) the average flux or count rate in any spectral band. The 
simplicity of the two shell approximation then allows to construct a large number 
of synthetic bursts to check if XRFs can be formed for some specific choice of 
the wind parameters. 



3. A statistical approach 

The redshift z and the four wind parameters r, E, T and k have been obtained 
in a statistical way for each synthetic burst. If long GRBs (and XRFs) are 
related to the explosive death of some special class of massive stars their birth 
rate is directly proportional to the star formation rate if) and their distribution in 
redshift can be deduced from tp(z) (Porciani and Madau, 2001). The distribution 
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Figure 1. Distribution of E p for the synthetic burst population: left 
(x = 1/2 and y = 5/2); right {x = y = 1/4). The dahed line represents 
the whole population and the full line the sub-group of bursts which 
can be detected by BATSE. 



of the observed duration tgo for long bursts is log-normal with a maximum at 
tgo ~ 20 s. We therefore also take a log- normal distribution for r (with a 
maximum at r max = 10 s) assuming an average burst redshift (z) ~ 1. The last 
three wind parameters E, T and re are very poorly constrained and we adopt for 
them uniform distributions: between 50 and 53 for Log-E, 100 and 500 for F, 
and 1 for Log re. 

We performed a first numerical experiment with 10 6 synthetic bursts to 
obtain their E p distribution. The results are shown in Fig.l for two choices of 
x and y: x = 1/2 and y = 5/2 i.e. synchrotron emission with standard equipar- 
tition assumptions and x = y = 1/4. This last case was already considered by 
Daigne and Mochkovitch (2003) who have shown that it leads to very good fits 
of the temporal and spectral properties of individual GRB pulses. The dashed 
line in Fig.l represents the distribution of E p for all 10 6 bursts while the full 
line corresponds to the sub-group of bursts which would have been detected by 
BATSE (a threshold of 0.2 ph.cm _2 .s _1 in the 50 - 300 keV energy band was 
assumed). It can be seen that for x = 1/2 and y = 5/2 the distribution of 
E p is much too wide while the agreement with the observations is excellent for 
x = y = 1/4. Figure 1 also shows that a large fraction of events with E p < 100 
keV is not detected by BATSE but may contribute to the population of XRFs 
seen by Beppo-SAX and HETE 2. 

To test this possibility we have generated (assuming x = y = 1/4) 1000 
events for which we compute both E p and the photon flux in the 50 - 300 keV 
band. The results are compared to the BATSE + Beppo-SAX data (Kippen 
et al, 2002) in Fig.2. Small dots represent BATSE bursts while large dots are 
events which are not detected by BATSE but have an X-ray flux (2 - 10 keV) 
larger than 1 ph.cm .s _1 . The similarity between the two diagrams clearly 
shows that XRFs can be produced in the context of the internal shock model. 
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Figure 2. E p - photon flux (50 - 300 keV) diagram: BATSE and 
Beppo-SAX data from Kippen et al (2002) (left) compared the result 
of our simulation (right). 



We then obtained the distributions of z, r, E, T and re for synthetic XRFs. 
Their redshift distribution is nearly identical to that of GRBs and our simulation 
therefore confirms that XRFs are not GRBs at very high redshifts. The duration 
and injected power distributions are also not very different between XRFs and 
GRBs and it finally appears that a reduction of the contrast re = r2/ri and an 
increase of the average Lorentz factor T are the most efficient ways to produce 
XRFs. The distribution of re has a sharp maximum at about 1.5 in XRFs while 
it steadily increases (and is negligible for re < 2) in GRBs. There are also five 
times more XRFs with f = 500 than with f = 100. 



4. Conclusion 



Our simulations have shown that XRFs can be produced by internal shocks if 
the Lorentz factor in the relativistic wind has variations of small amplitude only 
and a high average value. XRFs are not the result of a low injected power E but 
are weak and soft mainly due to a small contrast of the Lorentz factor which 
leads to a reduced efficiency of the dissipation in shocks. 
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